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SUMMARY 


It  Is  of  Interest  for  several  reasons,  from  both  the  mathe¬ 
matical  and  physical  points  of  view,  to  discuss  In  detail  what 
happens  to  the  various  categories  of  transport  equations  derived 
from  different  applications  of  invariant  Imbedding  as  the  velo¬ 
city  of  the  neutron  Is  allowed  to  become  arbitrarily  large  with 
a  corresponding  increase  in  the  probability  of  a  collision. 

This  idea  is  a  quite  natural  one  and  one  that  has  been  pur¬ 
sued  by  a  number  of  different  investigators  with  different  alms 
in  mind.  Diffusion  theory  classically  has  been  regarded  as  an 
approximation  to  the  more  rigorous  (but,  of  course,  not  completely 
rigorous)  transport  theory  under  the  assumption  of  high  velocity 
and  small  mean  free  path.  Purthermore,  passage  to  the  limit  in 
the  "telegrapher's  equation,"  a  linear  partial  differential 
equation  of  hyperbolic  type,  has  been  carried  out.  We  shall 
discuss  one  aspect  of  this  below. 

Our  principal  aim  here  is  to  study  the  limits  of  the  non¬ 
linear  functional  equations  obtained  from  the  transport  processes 
with  finite  velocity  as  the  velocity  Increases  without  bound. 

In  this  way,  we  obtain  corresponding  results  for  heat  or  diffu¬ 
sion  processes,  where  the  physical  picture  Is  not  as  clear. 

Having  obtained  the  equations  in  this  Indirect  and  complex 
fashion,  we  can  then  interpret  them  in  such  a  way  as  to  be  able 
to  derive  them  directly  by  invariant  Imbedding  techniques.  In 
all  cases,  the  equations  are  of  the  generalized  Riccatl  type 
which  we  recognize  as  characteristic  of  these  processes  of 
mathematical  physics. 
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INVARIANT  IMBEDDING  AND  NEUTRON  TRANSPORT  THEORY— V: 
DIFFUSION  AS  A  LIMITING  CASE 

Richard  Bellman 
Robert  Kalaba 
G.  Milton  Wing 


1 .  Introduction 

In  previous  papers  In  this  series,  we  have  Investigated  a 
variety  of  simple  models  of  transport  theory  by  means  of  the 
functional  equation  technique  of  Invariant  Imbedding.  A  number 
of  references,  together  with  some  discussion  of  more  complex 
versions  will  be  found  In  [l]  . 

Neutrons  are  mathematically  abstracted  to  be  point 
particles  with  finite  velocities.  Fission  and  scattering  are 
characterized  by  certain  probabilities  of  branching  and 
reversal  or  reorientation  of  direction  ("collision  cross- 
sections")  in  the  medium  within  which  the  process  is  occurring. 
In  the  great  proportion  of  cases  we  assume  no  neutron-neutron 
Interaction,  although  we  have  discussed  this  phenomenon  In 
one  paper,  [8]  . 

It  Is  of  interest  for  several  reasons,  from  both  the 
mathematical  and  physical  points  of  view,  to  discuss  In  detail 
what  happens  to  the  various  categories  of  transport  equations 
derived  from  different  applications  of  invariant  imbedding  as 
the  velocity  of  the  neutron  Is  allowed  to  become  arbitrarily 
large  with  a  corresponding  Increase  In  the  probability  of  a 
collision. 

This  idea  Is  a  quite  natural  one  and  one  that  has  been 
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pursued  by  a  number  of  different  investigators  with  different 
aims  in  mind.  Diffusion  theory  classically  has  been  regarded 
as  an  approximation  to  the  more  rigorous  (but,  of  course,  not 
completely  rigorous)  transport  theory  under  the  assumption  of 
high  velocity  and  small  mean  free  path,  [2]  .  Furthermore, 
passage  to  the  limit  in  the  "telegrapher's  equation,"  a  linear 
partial  differential  equation  of  hyperbolic  type,  has  been 
carried  out.  We  shall  discuss  one  aspect  of  this  below. 

From  another  direction,  the  discrete  random  walk  process 
yields  the  diffusion  equation  in  the  limit.  This  observation 
has  been  made  the  basis  for  a  considerable  amount  of  analytic 
and  computational  effort,  centering  about  the  theme  of  "Monte 
Carlo"  techniques. 

Cur  principal  aim  here  is  to  study  the  limits  of  the  non¬ 
linear  functional  equations  obtained  from  the  transport 
processes  with  finite  velocity  as  the  velocity  increases  with¬ 
out  bound.  In  this  way,  we  obtain  corresponding  results  for 
heat  or  diffusion  processes,  where  the  physical  picture  is  not 
as  clear.  Having  obtained  the  equations  in  this  indirect  and 
complex  fashion,  we  can  then  Interpret  them  in  such  a  way  as 
to  be  able  to  derive  them  directly  by  invariant  imbedding 
techniques.  In  all  cases,  the  equations  are  of  the  generalized 
Rlccatl  type  which  we  recognize  as  characteristic  of  these 
processes  of  mathematical  physics. 

At  the  present  time,  we  are  studying  the  question  of 
treating  Stefan-type  diffusion  problems  by  a  similar  passage 
to  the  limit  in  the  equations  derived  from  transport  processes 
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with  variable  boundaries.  This  is,  as  might  be  expected,  a 
complex  problem.  Some  initial  results  are  given  in  [9]  • 

Throughout  the  paper,  we  use  a  simple  generalization  of 
the  Idealized  one-dimensional  rod  process  treated  in  [6]  *  In 
the  following  section,  we  shall  obtain  some  new  equations  for 
the  flux  within  the  rod,  assuming  finite  velocities  initially. 
In  §3>  we  derive  Pick's  law  for  this  simple  process.  This  is 
important  for  our  purposes,  since  it  is  the  analysis  of  this 
result  which  suggests  the  combinations  of  functions  which 
should  be  used  in  the  limiting  case.  In  §4,  we  study  the 
limiting  form  of  the  internal  flux  as  the  velocity  becomes 
infinite,  and  in  §5  the  diffusion  process  giving  rise  to  the 
function  obtained  in  this  way  is  analyzed. 

We  then  turn  to  the  primary  objective  of  this  paper,  the 
passage  to  the  limit  of  the  nonlinear  integro-differential 
equation  obtained  for  the  reflected  flux  in  the  neutron  trans¬ 
port  case  by  means  of  the  technique  of  invariant  imbedding. 

In  the  concluding  section,  we  show  how  to  obtain  the  result  by 
direct  application  of  the  Imbedding  technique  to  the  diffusion 
process. 

Throughout  this  paper,  our  methods  are  largely  formal, 
since  we  are  principally  Interested  in  demonstrating  the 
applicability  of  invariance  principles.  The  existence  of 
relevant  limits  and  the  applicability  of  Laplace  transform 
methods  are  taken  for  granted  in  order  to  arrive  quickly  at 
the  desired  equations.  In  the  near  future  these  questions  will 
be  studied  in  a  rigorous  fashion. 

See  also  [3,4,5]. 
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For  tne  simple  mathematical  models  considered  hero,  there 
is  little  difficulty  in  carrying  out  thia  program.  Since, 
however,  we  know  that  the  passage  to  the  limit  involve a  a 
reduction  from  a  hyperbolic  partial  differential  equation  to  a 
parabolic  partial  differential  equation,  involving  inter  alia 

4 

a  redundancy  in  the  initial  conditions,  we  can  expect  some 
difficulties  in  the  general  case.  The  corresponding  study  for 
ordinary  differential  equations  when  a  limiting  value  of  a 
parameter  r?sults  in  a  drastic  change  in  the  order  of  the 
equation  is  of  sene  subtlety.  Particularly  interesting 
examples  of  equations  of  this  nature  occur  in  various  hydro- 
dynamical  investigations  where  viscosity  plays  the  role  of  th* 
parameter  which  approaches  zero;  cf.  Wa3ow,  [id]  . 

2.  The  Transport  Equation 

To  begin  our  work  it  is  necessary  to  write  down  some 
transport  equations  in  fairly  general  form.  While  3ome  of  them 
are  not  to  be  found  in  the  literature,  they  may  be  readily 
derived  by  the  methods  of  previous  papers . 

Consider  a  rod  of  material  which  transports  neutrons,  and 
let  the  neutrons  have  constant  velocity  c  (monoenergetic 
case).  The  usual  collision  processes  take  place  with  the 
probability  of  a  collision  in  a  length  A  of  the  rod  ta^en 
to  be  dA  +  o(6)  where  d  Is  a  constant.  On  the  average, 

2k  neutrons  emerge  from  a  collision  inside  the  rod,  k  going 
to  the  left  and  k  going  tc  the  right.  We  take  the  rod  to 
extend  from  0  to  x  (see  Figure  1),  and  designate  the 


u 


v~ 


coordinate  of  an  Interior  point  by  y . 
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Tc  initiate  the*  .:-r-cesa,  w;  auproae  to. at  there  if  a  tii.e 
dependent  .i.-.  e*,  ;(t  r.e;tia..^  Ir.e idem  to  the  icf*  pe: 


vj^iLouu  .1 


l  articled  at 


nd  none  at  t;.o  end  0.  Finally,  we  suppose  that 


We  write 


A 


c  r  x. 


v  0  i.  a  .  i  i  v  h  «  t  .  •  v  »  S*.  .  •  • 


(2.1 


f  *  ■  ♦  t  lin  j  v* .  -  /  ,  »  »  1  '  :*  r*  .  «'  p  -*lt  »  t'  **  '  j  ,1  •>  i 

\  J  f  ,  t  .  4  '  .  «  \l  W»  •  I  *•  *  »  h  .  d  t  »  O  I  W  j  •  V,  «.  M  .  tiw» 

*  ;t  n  O  1  £t  *  •  .  t  ♦"  !  **  I*  t  O  n(  *  *r  4  r^i  t  t 

puue  j»  u  t  w  X«  »•».  i*  e-. »  1  titv  v  Ailij  o  Vy  v  ■  ie< 

rii'it, 

(y,t  -  tne  average  nu:::ber  of  neutrena  per  second 


>  p  1  ',  ■■  *  •  •  ;.  •■  !  >r  .1 


i  --..I  >i  -  ■»  iv  .»  r 

Wk.  le*  Uiw  •  ^  *{^  C 


it.  * 


Using  tne  jr.etr.eaa  outlined  in  [a],  it  i3  easily  feund  that 


i  r.  o 

\  *-  •  *- 


*y 

*v 

<T7 


i 

c  .H 


+  — 


i  jv 

1  4- 

c/  w. 


-  v  iV  —  1  ,4  ♦  ^  \  , 

dVctl  4-  6  \  A.  -  1  V. 


V 


•t  -  1,  ^  \  4.  ,  *  -  ^  \  v  I  t  p>  C  , 


a(y,Ch  =•  v  (y  ,0 1  =0,  0  s  y  <  X. 


For  ac.ne  f  jr,.^;5  it  is  e  jivy  anient  to  talk  at  :-l  one  totfi  fi^.t 
from  tirj  C  tv  t.  *e  3nsix  tons,  l  slur.  *iy  use  capital  letters 
tv  inalcate  ^aantLtlea  integrated  eve:  tla.e.  Thao, 


r 


v2.3'i  U(y,t  =./)t  u(y,7.)di:, 

"0 


Q(t)  l(z)dz, 


etc.  It  Is  easy  to  see  that  tne  Integrated  quantities  satisfy 
equations  identical  to  '2.2  witn  the  low^r  case  .etters  belay 
replaced  by  capitals. 

Consider  n^w  the  case  in  whlcn  tne  source  consists  of  a 
single  "trigger"  neutron  at  t  =  0.  Taus ,  formally,  q(t;  = 
where  6  Is  the  Dirac  delta  function.  We  focus  attention  on 
the  particles  reflected  from  the  rod  at  x,  writing  w(x,t 
for  the  number  emergent  per  second,  and  W  ( x ,  t  for  the  total 

number  emergent  up  to  tine  t.  Clearly, 

(2.4  u(x,t ;  2  w^x,t  , 

U  { x  ,  t  '  s  W  ( x ,  t  . 


However,  It  is  well  it  regard  tne  x  in  the  arguments  cf  v 
and  W  as  referring  to  the  lengtn  of  tne  red  rather  thsn  to 
the  coordinates  of  the  end  point  of  the  r<-d.  With  this  rather 
subtle  distinction  in  mind  one  then  finds,  using  the  methods 


(2 


y  W  2  ?W 
?\  c  Jt 


dk  +  20  (k 


w  ( x  ,  o  W  '  x  ,  t 


z;a 


A  ;  X  ,  0  = 


w;o,t  =  o 


r  r. 


P-U:>. 
1  l-9-o9 


Notice  that  this  characterizes  tne  reflected  f^-x  in  a  fashion 
independent  uf  tne  interna,  fluxes. 

Pur  the  corresponding  case  in  which  there  is  a  source 
q(t),  we  write  F(x,t  for  the  reflected  flux.  Since  tne 
^  fundamental  physical  process  is  additive,  as  a  consequence  ol* 
our  tacit  sssumptior  that  tnere  are  no  interactions  between 
neutrons  passing  Ir.  opposite  iirectlons,  we  can  write 

(2.o)  P{x,t  • '  i<2  W(x,t  -  z'dz. 

C 

To  find  ar.  equation  satisfieu  by  F,  we  utilise  the 
Laplace  transform,  writing 

(2 . 7  P.  (x.s  «  Ax  e_stP(x,t  .dt, 

h  C/q 


with  a  consistently  similar  notation  fur  translurnis  of  ctner 
functions.  Then,  front  (2.~>  , 


(2 


26 


1  W.  C'ksW.  ~  , 

Li  Lt 


W.  ( 0,s  =  0 . 

Prom  (2  ,fc  ,  , 

(2.9  Pr  =  irw,  . 

L,  L> 

Hence  , 

'2-10  *  T  -L  +  2a<*  -  1PL  +  3kaFL-L' 

wnlcn  leads  back  to 


A- 
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(2.11)  *  f  p  OKQ(t)  f  2o(k  —  l)F(x.t) 


«•  ok  />t  w(x,z)F(x,t  -  z)dz. 


P(o,r )  •-  f(x,  ))  -  o. 

This  clearly  reduces  to  (2.8}  when  i(t)  »  8(t). 

fie  shall  derive  or.e  other  special,  case  of  (2. 11),  corres 
ponding  to  the  caae  when  i(tj  -  1..  *Jh'le  conceptually  this 
is  a  bit  harder  to  consider  than  the  single  trigger  neutron 
case.  It  has  the  mathematical  advantage  of  ;».v  '  ling  the 
S~ function.  For  this  type  of  source  we  write  the  integrated 
flux  as  r-:(x,t)  and  (2.11)  becomes 


(2.12)  ~  f  ~  -  okt  ~  2o(k  -  l)p(x,t: 


-  ok  ft  *  w(x,z)P.(x,t  —  z)dz. 


3at ,  from  (j.c  , 


(2.1?) 


P(x,t)  --  p r  w(x,t  -  z)dz 

if  r 


Then  we  easily  I'li.i 


(2.JU) 


f.  _1  o<t  ♦  ?a(>:  -  i  )  R(  x ,  t ) 


f  ok  ft "  r(x,z'r(x,t  —  z)dz, 

C/r, 


H(0,t)  *  R(x,?)  -  0,  r(x,t)  -  t~(x, t ) 

.  *  b 


5,  FIck  1  s  Law 
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If  we  subtract  tne  second  equation  of  (2. a  ircm  tne  first 
we  obtain 

(3-1,  ^(U4V  +J^(u~v  »-8(u-v. 

Per  large  c  ^ne  expects  tne  second  terrr.  on  tne  reft  to  be 
small.  Hence  we  formally  obtain  the  relation 

(1-2  +  V  -  -  d(u  -  V 

in  tne  "lln.it  of  large  velocity."  Equation  (/ .2  lc  ordinarily 
referred  tc  as  Pick 1 s  Law,  [2] ,  wnlcn  abates  that  tne  net  flux 
proportional  to  t..e  gradient  of  the  concentration  anu  Ir.  the 

up;.  ^tslle  direct  ton. 

4  .  Too  Ll.r  !_*.  leg  Co .  1  ’-ta  it  ted  fr  ‘2.2 

To  obtain  ra  1  ir.lna ry  results  we  take  Laplace  trans forme 

of  (2.2).  Thus,  using  the  notation  Introduced  In  (?./'), 
dU,  _ 

(4.1)  4  -  UL  =  o(’<  -  1  )u.  *  okvr , 

lv'  -  .  ,  ’ 

-  w  *  ~  ”1.  *  oicuL  ’  ~  l)V 

t*.  (■  ,F)  V  (x,s)  q.  (a). 

i-  u 


After  rather  extensive  but  rudimentary  calculations  we  get 

kcq.  (s  '  sinh  ?ly 


(4.2) 


ui.(y*8>  = 


A  cosh  Tlx  -  (■-  ♦  ( 1  —  k)o.sinh  ‘Axj 


1  (s}v/  cosh  7*y  t  (|  -r  (l  -  k)cr)slnh  ^y j 
v  (y,s  =•  — : - L - - - - - u. t 

>  - »-  O-  Or  j 


cosh  Tlx  -  +  (1  -  k)o)slnh  Tx 


with 


r-iSi? 

A  i  —  4 — 9 

10  ‘ 


(«.3) 


'  c 


x  r.Lt  --Ki. 


c  (1  -  2k) 


we  now  choose  k  -  i,  ?,  wh’ch  :ieins  physically  that  tii 

average  collision  ^Ives  rise  to  one  neutron.  This  choice 

eliminates  the  las"  tem  in  ( 4 .  j } .  Since  w«r  a- ek  1  dlffuiion 

type  equation,  we  let  c  — ♦  00,  l.e.,  the  velocity  become  infinite. 

Clearly  to  preserve  the  process  we  must  tnen  require  tnat  0  — *  uc 

In  such  a  way  tnat  11m  c/d  *  D,  a  constant.  Hence,  iin  A«=  '/s/D . 

(it  should  te  noted  that  a  somewhat  more  general  result 

could  have  teen  obtained  by  re^uirin^;,  ins.ea.  I,  -  1/2, 

2 

that  lim  a  (l  —  2k)  -  a.  by  so  doing  we  could  have 
accounted  for  oases  of  absorption  cr  ficoion.  Tc  Jo  this  here 
would  merely  coup.  J  oat*'  the  ensuing  calculations.) 

Bearing  ( 3 . 2’  in  mind,  we  set 


(4.4)  y  (y,$)  a(u.(y,3)  -  V  (y,s  , 

T  Xj  1/ 


i-im  J,  (y , 

o — 


Let  us  consistently  reserve  the  subscript  zero  to  refer  to  quan¬ 
tities  In  the  limit  .  s  *;  —>  a  . 
w'e  then  discover  that 


oqT  (  s  )  TV  ccoh  Ay 

^  |  Aw 

a — wd  ;  /.  coari  Ax  +  (-^  -*  ^)slnh  T.n 

r r_-f ", 

.  (S)  VI  cosh  (y»L  s 

s  1  nh  (  xJd  "^3  ) 


(4.5) 
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A  Classical  Diffusion  Problem 


..e  now  seel<  an  ordinary  diffusion  problem  which  gives  rise 
to  tne  limiting  expression  found  in  (^.5).  It  Is  reaully  veri¬ 
fied  that  if  ®(y,t)  is  Implicitly  determined  by  the  relations 

TT»  *  0,  «(x,t)  =  2q(tj,  G(y,C)  =  0, 


(3.1) 


>v?  3 1  ‘ 


then,  explicitly. 


(b.2>  er(y»3)  = 


2q.  ( s  '  t  i  nh  ( y  /d  ) 

Li 

3inh  (  X'/d  *  3  ) 


anr 


(5. a) 


dy 


2q.  (s  )VD  Js 


l  n”1 «  SS^LlIrIIi 


rrr" 

sinn(xVr  s ) 


.*e  may  summar* /e  our  results  thus  far  as  follows: 

If  we  consider  the  transport  problem  formulated  In  (f. .2) 

In  the  ilii.'tlng  care  where  c  — ►  oo  ,  c/a  — *  D,  with  k  =•  i/2, 
then  the  ^rob-iem  Is  formally  equivalent  to  the  classical  dif¬ 
fusion  prcolen  (b*l  '•  The  quantity  l*m  (u(y,t)  +  v(y,t') 

j—*oo 

may  be  Identified  with  w(y,o'j,  whl le  11m  d(u(y,t)  —  v(y,t)) 

. .  .  0_^00 

corresponds  tc  —  3 e  '$ y . 

It  is  possible  tj  Identify  €»(y,t  with  the  total  neutron 
flux  (see  [  2  ]}  although  the  diffusion  may  refer  a3  well  to 
heat  or  material  cor.  .ent  ration.  The  fact  that  a  sourc®  c  f 
2q(t)  is  requires  in  trie  problem  (5.  i  )  :ruy  be  rather  puzzling 
unt’i  .Tic  notes  r  >r  ( k .  s.  hat,  forr..  i  iy ,  octh  u(x,t)  and 
v(x,  c  )  apprvac.u  (  t)  vs  c  — >  cc  . 


acted  F.ux 


P-  • 

n-o — o 

12 


p  13  j  f' 

^  •  4t«V.  *  l  4  . 


o-et  us  now  turn  t  liquation  ic'.n  :  ana  try  to  carry  <.ut 
tr.e  saoie  type  ..f  pasca^o  tv.  t ne  liii.it.  It  Is  clear  that  we 
si.ust  begin  by  Investigating,  tr.e  qu-ritity 


i.i 


H%x,r  =  .<R(x,t  -  Q't  r  . 


which 

reduces 

to 

cO 

R(x,t  ' 

—  t ;  , 

since  q(t 

=  i .  Triua 

(  .2) 

M  X , 

t.  j 

a 

uf  x  r.  ’ 

0 

♦  t. 

r(x. 

t } 

h 

(AlU 

0 

1. 

3ubst 

1  tut 

i  n.-  * 

hesr  1 

r.  ( i .  - 

)  ' 

i  -  * 

find 

(6.3) 

: 

0  )T 

• 

(ii 

?  '  c 

~  1) 

<yt 

/  H 
°(o 

- 

+  .2  /? 

t  + 

h(x,t  -  2) 

:v/c 

(1 

0 

Pro  7i 

this 

we  re -id 

-  •  *, 

ti** 

/  ii ' 

V  •  4  l 

3.3  . 

0  X 

C 

(-0 

1  I)  = 

1  /?* 
2< 

',(v,2h(x.t 

-  z  ;  d  l , 

M{x, 

? 

-  0 

.  H(. 

t  —  . 

#  /  “ 

-  cc . 

t< 

•)  the 

ill 

ml  t 

] 

vS  ,  1  n 

pe 

^efc  (at  i‘is 

t  *'o  rma  i )  y  ) 

ib.5) 


Fx 


..r-i 

+  <iv  * 


H0(x,0)  -  0,  Hn(0,t 


(x#2}h0(x,t  -  Z) dz, 

i 

=  r  «;  1 


~K  ' 
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where  H,(x, t )  iir?.  H(x,t),  *tc.,  3 is  agreed.  That.  (6.5)  is 

c— > oo 

the  correct  limiting  form  may  be  established  by  a  Laplace 


trans  ’  rv.  argument  similar  to  that  of  the  last  section.  We 


unit  the  let 3 12  . 


It  's  of'  some  interest  to  evaluate  H„.  This  may  be  done 


by  solving  (’.5),  c:  course.  Hcwever,  it  is  easier  for  uu  to 


no t  e  ’’  rc.  r  ( k .  5 )  that 


(  -6; 


h.  (x,s)  *  -  2q  ( .s 


s’.nh  ( x  Jd  \ ) 


— —  co  tii  (x/r  ‘si. 


since  q(t)  -  i.  We  find 


2  2 

x  n 


(•  h,  (x,t; 


;  -t  ^  e 


f 


/  i  t 

~  x * 

JX 


w.ner*  -  .  theta  ’met 'on. 


The  i-i<tlogu**  (  .5)  '•••y  be  ie  rive  a  easily  for  the  case 


' n  which  ther-  Is  an  irb*  t  nary  sconce  : ( t ; .  The  reeult  1» 


(6.0)  *  21  J  -  4(//?‘  M*»:  )Jc(x»t  ”  z)dz. 


J,.{x,0)  »-  >,  J^(C,t) 


•  1  ’  rr.  c  )F(x,  *  )  -  W.(t)/# 


where  F 


Is  as  in  (2.6). 


V 


now  read*  j  s*»c  the  roii owing  rt  jit  . 

T. .  consider  the  r-..-sp  jgroc  . _ r  .  |ted  in  ( 2.2) 

lit  :ne  . ‘mlt'ng,  case  then  the  :u.-.v  1  •- ./  J  ( x , '  1_ _  formal xy 

e  jy.l  va.ent  to  the  Quart I ty  ~  /'^  0(y,3)ds!  where  w 

y( 0  1  y  =  x 

1 3  define.:  ,y  ( „  .  i  } .  Purl  her  J  ,  a  w  1 1 3 ;  len  ( \  .  b )  with 
h ,-,(*#*  ;  J  -er.  (.;■  (b.5). 


7.  A  aired  Invariant  Lmte-da  lap  Aj-yroacn  In  Diffusion  Theory 

The  aquations  thus  far  obtained  tie  not  i.cn,  though  oar 
approach  to  them  may  be  somewhat  novel.  To  conclude  our  work 
here  <ve  si.al ;  present  a  nethol  '  ->bta’  nJng  (t  .8)  i\,  Ir.vu riant 
1  :noed  iir»v*  ques  without  ve:  * :  ring  outr’d-  the  confine  i  of 

or-  i  tar^  r 1  "  h  «  >r  .idry.  Th”  n-  3  3 e err l 1  ed  holds  promise 

o“  being  a.  .  I  D  .  ir.  at  .1  vore  11  .ab  )  31  f  "us!  on  t  r>o- 

ceas->-  than  that  ies'orlbed  here,  trd,  In  particular,  may 
e-  .t  /aa’.  iy  yield  new  born.. latl'ris  c.  te-  an-  ype  pr  lems. 

Tc<  be  c<  nai stent  n  c  r  viewpoint,  we  now  thin*  or  ,t) 


as  the  density  c  rie  ut  rone  at  y  at  tire  t.  Then  the  net 
neutron  current  ien?Ii  c ( y , t )  Id  prcvicieu  by  Flck's  law.  In 
the  ordinary  diffusion  approach  [2  j. 


,7.1  J 


O  ; 


The  conservation  cl*  part 
production.'  require-  n 


-  > 


C  (  3  ,  t 


A: 


1  ;  >8  (s'r.ce 

ar.j  Interval 


) 


J  /;b 

M 


there  is  r.  internal 
(a,b)  of  the  rod 

p(y » v  'd^  - 


( 


II 


.  .c?  ’  .  s  a  r.'i  l 

*  . 

or  t h n 

*»t  current  en.e 

rg'rg  “rom  our  roa 

.eugth  x. 

(  X  ‘ 

V  A  t 

Hf  e 

,  a  ga  l  n ,  .vhlle 

it  is  t rue  tha t 

r.  (  X  ,  - 

.%«• 

choose  • 

<>  ■'egard  the 

x  'n  the  '*inc  t  ion 

,  ^  •  referring 

t  -  t 

he  ler.^tn  0  ‘  the  rod. 

Thus  x(x  +  L,t' 

Is  the  ne;  rre 

nt  e 

’oerger.  t 

•*r>:r  ■  rod  or’ 

length  x  +  £>, 

f  u  rc  e  J  .  { ■ 

r 

V 

X  >  i  # 

orher  1 nt  ■  '3i 

and  Loundary 

:or.  i !  f  rt5  k  '  t  .g 

an 

e 

At?  iv«w  try 

t"  e 

•  re 

•<{>  -  .,,1,  IrVterrns  0!'  k(  x,  t v . 

A,f lying  ( 

'  t  h 

r*o  1  0  f 

iength  x  *  fi 

tvo  *'ind 

(..  1  4*  * 

9 

-  U. 

1  rfX+L 
,*  T  C/ 

X 

J»(y,h  )dy  , 

r,  ’  *1  teg  rat  ‘  n. ; 

ve  r 

*  *  t*  ^ 

-  • 

(  .*  r‘(x 

■  t 

C(x,t 

-  <*x  of 

^  X 

if  /  t  Jy . 

..e  r: '  .  : 

exj 

:•«?  >*  . 

and  d.  To  ind 

'  (  x  ,  ’ 

tha  - 

•»e  "  a  ve 

hi  0  ■  a’  disregarded  the  part  of 

t.t*  r  :  :‘r/; 

1 0 

x .  V 

Ir.e  :oi  U:u’ >y 

conditions  Impose J 

iJ  ;  .2;  .  .  ’ '"Hr 

■  .<  n  . 
%•  • 

»?  r.L'rt 

*hu  (x*t} 

is  r ereiy  the 

"u:  rent  o  re  r  •  .e  ;  the  .r  e  p(x, t }  imposed.  Let  us 
;  ’  steady  source  o;  unit  strength  produces  a 

r  rrvrit  out  {.'(it,:. ) .  The.*  a  source  p(x0*  ) 

.  produce  a:  integrated  current 

(  o!  C(x,  ^0'  *(x,t-  z(p(x,z)dz. 

(This  4  •  'us*  .  *avia:*«e  1 ' 5  rrinciyle  [  7  | ) . 


As  y~t  we  have  not  used  (7.1).  From  It  we  rind 


F -j  S  ;•; 


(7.6)  ^(x,t)  =  ±  k(*  *  A,t)  +  2q(t)  +  o(A). 


Substituting  (To)  and  (7.0)  In  (7. 4),  we  ootaln 


(7.7) 


:'.(x  *  i,t)  p(x,t  -  x)|^k(x  b, 2)  >  2q(z)jdz 

--  2Aq(t)  +  o(A). 


But ,  ti.v  Duhame  1  r  i  nr  ’  r  1  e  , 


4* 


(x,t  -  2)1(7: dz  K(x,t; 


Thus 


(7.9)  |tt  ■»  2-i(t)  p(**fc  -  z)k(x,z)dz. 


This  agrees  with  ( o . H ,  upon  identifying  k  with  DJ  and 

n 

i'  w’th  £  h„,  the  :'actor  1/2  occurring  because  p  Is  the 
current  due  to  a  unit  source,  while  hr  4s  obtained  1‘rom  a 
source  of  strength  ?. 

It  is  clear  that 

(7.10)  K(x.O)  *  0. 

T'*'  find  K(0#t)  w*  not-3  from  (7.r)  that 

«(0ft)  -  0  *  A  k(',t)  +  2q(t)  *  0 (L), 


30  that  for  *(t)  >  C, 
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